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ABSTRACT 


A one-dimensional  model  of  the  tropical  atmosphere 
consisting  of  conservation  equations  for  heat,  moisture, 
momentum  and  mass  is  constructed  and  used  to  investigate  the 
effect  of  sea  surface  temperature  perturbations  on  the  local 
thermodynamic  characteristics  of  the  marine  boundary  layer 
in  the  summertime  Western  Atlantic  trade  wind  regime. 

The  height  of  the  trade  wind  inversion  and  the 
thickness  of  the  atmospheric  mixed  layer  are  found  to  be 
little  affected  by  changes  of  sea  surface  temperature  in  the 
range  that  would  be  produced  by  Ocean  Thermal  Power  Plants. 
The  surface  fluxes  of  long-wave  radiation,  short-wave  radia- 
tion, and  sensible  heat  are  practically  unaffected  by  these 
small  changes  in  the  sea  surface  temperature  or  by  changes 
in  the  intensity  of  the  cumulus  convection.  The  surface 
latent  heat  flux  (taken  positive  downward),  however,  varies 
linearly  with  the  sea  surface  temperature^ wi%h-e^-&4ope  of 
ifemit — 105r-cffT  cin~^  day“^~  °K^ *<£and~Ts--strong  1 y coupled  to 
the  intensity  of  the  cumulus  convection. 

If  the  maximum  tolerable  change  in  the  surface 
evaporation  rate  is  taken  to  be  at  -5%,  the  maximum 
allowable  OTPP-generated  sea  surface  temperature  perturba- 
tion is  about  -0.  2^C. 


Section  1 


INTRODUCTION 


The  operation  of  one  or  more  Ocean  Thermal  Power 
Plants  (OTPP's)  in  a particular  region  will  cool  the  upper 
layer  of  the  ocean  and  reduce  the  sea  surface  temperature 
(SST)  by,  perhaps,  as  much  as  1.5°C  (Piacsek,  et  al.  , 1976). 
This  will  alter  the  net  surface  heat  flux,  EQ,  the  local 
thermodynamic  structure  of  the  atmosphere,  and  possibly 
the  climate  of  the  region  downwind. 

It  is  important  to  assess  these  environmental 
impacts  before  large  numbers  of  OTPP's  are  put  into 
operation.  This  point  is  underscored  by  the  work  of 
Shukla  (1975)  who  performed  experiments  with  a general 
circulation  model  and  found  that  an  average  SST  anomaly  of 
-1.5°C  in  the  Arabian  Sea  (over  a region  roughly  the  size 
of  the  Gulf  of  Mexico)  reduced  the  model  Monsoon  rainfall 
in  India  by  about  50^c.  This  relationship  between  SST  and 
rainfall  was  qualitatively  supported  by  Shukla  and  Misra 
(1977)  who  analyzed  60  years  of  SST,  wind  speed,  and  rainfall 
data  for  that  region.  Such  evidence  linking  rainfall  to  SST 
is  disturbing  since  the  Gulf  of  Mexico,  which  is  the  source 
region  for  most  of  the  summertime  rainfall  in  the  mid-western 
United  States,  is  also  a proposed  OTPP  site. 

In  addition,  Wendland  (1977)  analyzed  nine  years 

of  SST  and  tropical  storm  frequency  data  for  the  Atlantic 

and  found  that  decreases  in  SST  of  a few  tenths  of  a degree 

6 2 

Celsius  over  regions  with  areal  extent  of  order  10  km  , 
were  associated  with  a significant  decrease  in  the  number  of 
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monthly  storms  and  the  length  of  the  storm  season. 
Wendland's  findings  imply  that  large  numbers  of  OTPP's 
operating  in  the  Caribbean  may  decrease  the  tropical  storm 
activity  there.  Since  a significant  percentage  of  the 
annual  rainfall  in  the  southeastern  United  States  is  due 
to  tropical  storms,  it  is  not  clear  whether  this  would  be 
desirable. 


To  ensure  that  OTPP's  are  deployed  in  such  a 
manner  that  the  environmental  changes  they  cause  fall 
within  an  acceptable  range,  it  is  necessary  to  determine 
the  change  in  the  net  surface  heat  flux  with  respect  to 
sea  surface  temperature,  dEo/dTQ.  Piacsek  et  al.  (1976) 
and  Bathen  et  al.  (1976)  used  bulk  aerodynamic  formulation 
and  local  meteorological  data  to  estimate  dEQ/dTo  at 
proposed  OTPP  sites  near  Puerto  Rico  and  Hawaii,  respectively. 
With  this  type  of  approach,  however,  it  is  necessary  to  make 
arbitrary  assumptions  about  the  atmosphere's  thermodynamic 
adjustment  to  the  SST  perturbation.  Another  possible  way 
of  estimating  dEQ/dTo  is  to  compare  direct  measurements  of 
the  surface  fluxes  at  a particular  location  at  another  time 
of  the  year  to  those  made  at  the  same  location  at  a different 
time  of  the  year  when  the  SST  is  different  (summer  and  winter, 
for  example).  However,  the  large-scale  flow  patterns  of 
the  atmosphere's  general  circulation  system  also  vary  with 
the  season  and  this  would  alter  the  seasonally-averaged 
surface  fluxes  whether  the  SST  changes  or  not.  Consequently, 
there  would  be  no  way  to  identify  how  much  of  the  observed 
change  in  EQ  was  due  purely  to  the  change  in  SST.  A more 
reliable  approach  is  to  explicitly  model  the  atmospheric 
response  to  changes  in  SST. 
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A rather  sophisticated  model  is  necessary  for 
this  type  of  study  since  it  must  be  capable  of  reproducing 
the  observed  average  structure  of  the  tropical  marine 
boundary  layer.  This  structure  has  been  detailed  by 
Garstang  and  Betts  (1974)  for  the  undisturbed  (i.e.,  fair- 
weather)  trade  wind  regime  and  consists  essentially  of 
(see  Figure  1): 


• A surface  or  "constant  flux"  layer 
extending  from  the  sea  surface  to 
typically  50  m height  in  which 
potential  temperature  and  water 
vapor  mixing  ratio  decrease  with 
altitude ; 

• A subcloud  mixed  layer  extending 
from  the  top  of  the  surface  layer 
to  roughly  500  m height  in  which 
potential  temperature  and  water 
vapor  mixing  ratio  are  approximately 
constant  with  altitude; 

• A thin  transition  layer  at  the  top 
of  the  mixed  layer  in  which  poten- 
tial temperature  increases  and 
water  vapor  mixing  ratio  decreases 
sharply  with  altitude;  and 

• A conditionally  unstable  cloud  layer 
extending  from  the  top  of  the  tran- 
sition layer  to  about  1-2  km  height. 


The  marine  boundary  layer  is  topped  by  the  trade  inversion 
in  which  potential  temperature  increases  and  water  vapor 
mixing  ratio  decreases  rapidly  with  altitude.  Below  the 
transition  layer  the  vertical  eddy  fluxes  of  heat  and 
moisture  are  due  mainly  to  dry  convection  and  mechanical 
mixing,  while  above  it  these  fluxes  are  accomplished 
primarily  by  cumulus  convection.  Thus,  both  dry  convection 


and  cumulus  convection  provide  important  coupling  mechanisms 
between  the  sea  surface  and  the  upper  part  of  the  marine 
boundary  layer. 

Recently  a number  of  diagnostic  studies  have  been 
carried  out  which  clarify  the  role  of  cumulus  convection  in 
maintaining  the  thermodynamic  structure  of  the  tropical 
atmosphere.  In  these  studies  a model  of  the  cumulus 
ensemble  is  used  along  with  observed  large-scale  budgets 
of  heat  and  moisture  to  deduce  the  characteristics  of  the 
clouds  and  their  effect  on  the  environment.  This  method 
was  first  applied  to  the  undisturbed  trade  wind  regime  by 
Nitta  (1975),  who  used  a model  fundamentally  similar  to 
that  of  Arakawa  (1971)  to  analyze  the  BOMEX  Phase  III  data. 
He  found  that  most  of  the  clouds  were  confined  to  below 
the  trade  inversion  and  that  the  large  observed  heat  sink 
and  moisture  source  in  and  near  the  base  of  the  inversion 
was  due  primarily  to  the  cooling  and  moistening  effects  of 
the  detrainment  of  liquid  water  and  saturated  air  from  the 
clouds.  These  conclusions  were  also  indicated  by  the  work 
of  Soong  and  Oqura  (1976)  who  analyzed  the  same  data  set 
with  a model  consisting  of  six  cloud  types,  each  of  which 
were  simulated  with  a two-dimensional  axisymmetric  cloud 
model . 


Representing  the  effects  of  cumulus  clouds  in  a 
numerical  prediction  model  is  a difficult  task.  One  way 
to  accomplish  this  is  to  use  a very  high  resolution  three- 
dimensional  model  which  explicitly  resolves  each  individual 
member  of  the  cumulus  ensemble  and  simulates  its  full  life 
cycle.  This  approach  was  taken  by  Sonvneria  (1975)  in  a 
remarkable  paper.  However,  while  this  technique  gives 
great  insight  into  the  relative  importance  of  various 
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physical  processes  and  clarifies  the  relationship  between 
the  cloud  field  and  the  subcloud  layer  turbulence  field, 
it  is  practical  only  for  very  limited  studies  because  of 
the  tremendous  computer  resources  required. 

A more  common  approach  to  the  problem  is  to 
parameterize  the  effects  of  cumulus  clouds  in  terms  of 
the  model-predicted  large-scale  fields.  This,  of  course, 
is  intrinsically  a more  difficult  task  than  the  diagnostic 
determination  of  cumulus  ensemble  characteristics  since  it 
required  prediction  of  both  the  characteristics  of  the 
ensemble  and  the  time  evolution  of  the  large-scale  fields. 
Cumulus  parameterization  has  been  one  of  the  most  important 
problems  considered  in  tropical  meteorology  and  has  under- 
gone rapid  advances  in  recent  years.  Early  attempts  at 
cumulus  parameterization  consisted  of  simple  convective 
adjustment  schemes  which  conserved  energy  but  did  not 
explicitly  take  into  account  important  processes  such  as 
compensating  subsidence  in  the  environment,  entrainment 
and  detrainment,  and  the  re-evaporation  of  cloud  droplets. 
The  landmark  paper  of  Arakawa  and  Schubert  (1974) 
incorporated  these  processes  into  a spectral  model  of  a 
cumulus  ensemble  and,  for  the  first  time,  provided  a 
closed  and  theoretically  rigorous  parameterization  of  the 
mutual  interaction  of  a cumulus  ensemble  with  the  large- 
scale  environment. 

Modeling  of  the  tropical  subcloud  mixed  layer 
has  also  progressed  rapidly  in  the  past  few  years. 

Following  the  earlier  work  of  Ball  (1960)  and  Lilly  (1968), 
Betts  (1973)  proposed  a model  which  showed  how  the  thermo- 
dynamic characteristics  and  the  thickness  of  the  mixed 
layer  are  controlled  by  large-scale  mass  convergence, 
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cloud-induced  subsidence  in  the  environment,  and  the 
radiative  and  surface  fluxes.  Ogura  and  Cho  (1974) 
included  the  cloud-induced  vertical  eddy  fluxes  of  heat 
and  moisture  at  the  top  of  the  mixed  layer  into  a model 
similar  to  that  of  the  earlier  studies.  In  all  of  these 
models,  the  transition  layer  is  treated  as  a discontinuity 
in  the  thermodynamic  variables  at  the  top  of  the  mixed 
layer.  These  studies  are  important  in  the  present  context 
because  they  indicate  that  the  steady-state  surface  fluxes 
are  determined  by  complex  interactions  involving  the  entire 
marine  boundary  layer. 

Pond  (1971,1972)  reviewed  surface  layer  models 
useful  for  predicting  the  surface  fluxes  in  terms  of  mean 
field  quantities.  Since  then  a number  of  experiments  have 
been  performed  which  have  increased  our  confidence  in  the 
empirical  forms  of  the  Monin-Obukhov  similarity  functions 
which  arise  in  this  type  of  treatment.  Busch  (1977)  has 
given  a comprehensive  review  of  recent  developments  in 
this  area. 


A reliable  determination  of  the  atmospheric 
response  to  SST  perturbations  requires  a comprehensive 
model  which  can  realistically  simulate  the  complex  physical 
processes  responsible  for  maintaining  the  structure  of  the 
marine  boundary  layer.  Such  a model  has  been  constructed. 
The  remainder  of  this  paper  will  describe  it,  show  some 
preliminary  results,  and  indicate  plans  for  future  work. 
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Section  2 

DESCRIPTION  OF  THE  MODEL 

FUNDAMENTAL  EQUATIONS 


Consider  a region  above  the  tropical  ocean  large 
enough  to  contain  an  ensemble  of  cumulus  clouds  but  small 
enough  to  be  regarded  as  a fraction  of  a large-scale  system. 
Taking  pressure  as  the  vertical  coordinate,  the  equations 
for  conservation  of  heat,  moisture,  momentum  and  mass 
spatially  averaged  on  constant  pressure  surfaces  in  the 
region  can  be  written  (Yanai,  Esbensen,  and  Chu  (1973)) 


y*Vs  + 
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9P 
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Lv(c  - 

— . a — » — » 
e>  - “ 
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e - 

a - 
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f - 

3P 

0 
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where  s is  the  dry  static  energy  given  by 


s ■ CpT  + gz 
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q the  mixing  ratio  of  water  vapor  (mass  of  water  vapor  per 
mass  of  dry  air),  u the  east-west  component  of  wind  speed, 
v the  north-south  component  of  wind  speed,  j the  horizontal 

wind  velocity  vector,  p the  pressure,  u the  vertical  p- 
velocity  u>  « ^ , QR  the  net  radiational  heating  rate,  c 
the  rate  of  condensation  per  unit  mass  of  air,  and  e the 
rate  of  evaporation  per  unit  mass  of  air.>  Other  symbols 
are  defined  in  Appendix  A and  all  of  the  notation  is  standard. 
Spatial  averages  at  constant  p taken  across  the  entire  model 
domain  are  denoted  by  (-) , primes  indicate  departures  from 
these  averages,  gradients  taken  at  constant  pressure  are 
denoted  by  (V),  and  we  assume  that  fluctuations  in  the 
horizontal  components  of  wind  are  uncorrelated  with  s',  q', 
u ' , an  d v ' . 


We  also  assume  that  the  fraction  of  the  model 
domain  covered  by  active  cumulus  clouds,  o,  is  small.  Thus, 
since  s and  q inside  cumulus  clouds  over  the  tropical  ocean 
tre  not  observed  to  depart  substantially  from  s and  q of  the 
environment,  we  have 


s - s 

q ~ q , (7) 

where  (~)  denotes  horizontal  averages  taken  across  the  clear 
region  between  the  clouds.  Furthermore,  from  mass  continuity, 
we  have 

u)“iU_+  (1  - o)u>  ~ w + w , (8) 

o o 
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with  wo  given  by 


i -i  °i 


where  is  the  average  vertical  p-velocity  at  a particular 

level  inside  the  i-th  cloud,  oi  is  the  fraction  of  the  model 

domain  covered  by  the  i-th  cloud,  and  N is  the  total  number 

of  clouds.  Note  that  uL  is  observed  to  be  large  relative  to 

o 

w.  Thus,  in  general, 


b)  yt  U> 


In  addition,  we  have  the  hydrostatic  equation 


- Pg 


the  equation  of  state  for  moist  air 


£ = RT( 1 + 0.61q) 


(ID 


and  the  equation  for  potential  temperature 


6 - T 


W 


(12) 


where  p is  the  density  of  the  moist  air,  g the  acceleration 
of  gravity,  z the  altitude,  R the  gas  constant  for  dry  air, 
pQ  the  surface  pressure,  T the  temperature,  and  cp  the 
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specific  heat  at  constant  pressure  for  dry  air.  Note  that 

-2 

since  q is  at  roost  0(10  ),  we  assume  that  the  specific  heat 

of  moist  air  is  given  by  the  specific  heat  of  dry  air,  which 
is  common  practice.  This  along  with  the  assumed  hydrostatic 
balance  of  the  environment  gives  from  (6)  and  (12) 


= c fa 

3z  Cp^] 


V 


R/cp  36 

Tz 


(13) 


Thus,  a layer  of  constant  potential  temperature  is  also  a 
layer  of  constant  dry  static  energy. 


The  following  six  subsections  give  a comprehensive 
discussion  of  various  simplifying  assumptions  applied  to  the 
basic  model  equations,  the  parameterizations  used,  and  the 
numerical  techniques  employed.  The  reader  who  wishes  only  a 
summary  of  the  resulting  model  formulation  may  skip  to 
Section  2.8. 
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THE  SUBCLOUD  MIXED  LAYER 


2.2.1 


Thermodynamic  Equations 


Consistent  with  the  observations  of  Bunker  et  al. 

(1949),  Malkus  (1958),  and  many  others,  we  assume  that  a 

mixed  layer  exists,  extending  from  the  top  of  the  surface 

layer  to  a level  denoted  by  pb>  in  which  s and  q are  constant 

with  height.  We  denote  the  pressure  at  the  top  of  the 

surface  layer  by  pn  , the  pressure  at  cloud  base  by  p , and 

restrict  our  attention  to  the  typical  trade  wind  situation 

whre  p„  > p.  > p„  in  all  that  follows, 
a b c 


We  will  model  the  subcloud  mixed  layer  with  the 
formulation  of  Arakawa  and  Schubert  (1974).  In  this  type 
of  treatment  the  transition  layer  is  modeled  as  a discontin- 
uity in  r.  and  q at  pb-  This  is  illustrated  in  Figure  1. 

Defining  s_  and  q as  the  horizontally-averaged 
m m 

mixed  layer  values  of  dry  static  energy  and  water  vapor 
mixing  ratio,  and  assuming  no  condensation  or  evaporation 
of  cloud  droplets  or  raindrops  occurs  in  the  mixed  layer, 

(1)  and  (2)  can  be  integrated  over  the  mixed  layer  depth  to 
yield 

3s  , v [CTT’L  " 0^>. I v 

n®  + (r’O  - J j^r L * («r)u  ■ 


(14) 
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y.vq  . 1<»«  >b-  <“«  >a. 

3t  x M Pb  - pa 


(15) 


where 


1 /D 

< >M  S pT^T  j < > dp 

pb  pa  p. 


The  terms  (w's')b  and  (w  q )b  are  the  turbulent  fluxes  of  s 
and  q at  the  base  of  the  transition  layer.  They  arise 
because  the  turbulence  at  the  top  of  the  mixed  layer  due  to 
dry  convection  entrains  the  generally  warmer  and  dryer  air 
above.  As  mentioned  in  Section  1,  it  is  observed  that  the 
eddy  fluxes  due  to  dry  convection  are  confined  to  below  the 
top  of  the  transition  layer.  Thus  the  entrainment  fluxes  go 
to  zero  at  the  top  of  the  infinitely  thin  transition  layer 
which  is  denoted  by  p^. 

We  now  define 


As  = s(pE)  - sm  , 


Aq  i q(PE)  - Qm 


(16) 


and  consider  the  heat  and  moisture  budgets  in  the  transition 
layer.  The  total  downward  fluxes  of  s and  q through  the  top 
of  the  transition  layer  in  the  region  between  the  clouds  are 
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> 


(s  + As) 
m 


3pb  - — 

TT  + Vvpb  + “b 


(qro  + Aq) 


Jb  . ^ ^ 


+ ^b*Vpb  + u>b 


(17) 


while  the  downward  fluxes  of  s and  q through  the  bottom  of 
the  transition  layer  in  the  cloud-free  region  are 


- (al6  * *(>•%)  * 2b  + <"  s >i 

- + V’Pb)  + ^b  + <“^>1 


(18) 


where  wb  is  the  average  vertical  p-velocity  in  the  environ- 
ment between  the  clouds  at  pb>  Combining  (17)  and  (18)  and 
requiring  the  fluxes  of  s and  q to  be  continuous  across  the 
transition  layer  yields 


>b  . = , <=T>b 


at-  * *b-7»b  ' % - 


3Pb  - — 

+ V^-VP,, 

3 1 -*-b  b 


wb  - 


(^TTr)b 

Aq 


(19) 


which  imply 


(wT*q"T  )b 


t§  <^>b 


(20) 
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The  horizontally-averaged  virtual  dry  static  energy 
is  given  to  excellent  approximation  by 

s v = s + 0.61  cp  T q , (21) 

where  the  term  arising  from  the  correlation  of  T'  and  q' 
has  been  neglected.  The  virtual  dry  static  energy  accounts 
for  the  effect  of  moisture  on  the  buoyancy  of  an  air  parcel. 
Following  Lilly  (1968),  Betts  (1973),  and  others,  we  assume 
that  the  turbulent  flux  of  virtual  s at  the  top  of  the  mixed 
layer  can  be  written 


(=^)b  = - k(^)0  • <22) 

where  / w',s.M  is  the  surface  flux  of  s , given  to  excellent 
\ v/o  v 

approximation  by 


(ST^s l ) = (uTri"T)o  + 0.61  cp  T(^T^)o  , (23) 


and  k is  a constant.  Equation  (22)  arises  from  considera- 
tion of  the  turbulent  energy  balance  in  the  mixed  layer 
which  requires  k to  be  in  the  range  of  0 to  1.  Deardorff 
et  al . (1969)  suggested  k ■ 0.10  from  laboratory  experi- 
ments, while  Betts  (1973)  suggested  k * 0.25. 


Combination  of  (20)  through  (23)  yields 

<=^>b  - - ^ (=^)0 

<sV)b  - - (“^)o 


1 
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where 


Asv  « As  + 0.61  cp  T Aq  . (25) 

In  general,  Aq  is  negative  while  As  and  Asy  are  positive. 
Note  that  Asv  > 0 is  required  for  dry  convective  stability 
at  the  top  of  the  mixed  layer. 

Equations  (8),  (14),  (15),  (19),  and  (24)  can  be 
combined  to  yield 


W5  - - - 57b;  [<“*>.♦  kH;(“  sv)J  + («*)„ 


(26) 


3q, 

W 


m 


- * 5-r-^  [<“  Q )a  + kif;(“  sv) 


(27) 


3pt 

TT 


- vVpb  * % - K)b  * zi;(“  sv)( 


(28) 


which  are  the  equations  which  model  the  thermodynamics  and 
thickness  of  the  subcloud  mixed  layer. 
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2.2.2  Momentum  Equations 

We  will  assume  that  the  clouds  have  no  direct 
effect  on  the  vertical  eddy  momentum  flux  in  the  mixed 
layer.  In  addition,  we  will  assume  that  the  Rossby  number 
of  the  horizontally  averaged  flow  is  <<1,  allowing  us  to 
write 

H - - = H ♦ *(*  - *,)  - ^ • <29» 

If  ■ - “ If  + f(“g  - • <29b) 

which  come  from  (3)  and  (4)  with  the  horizontal  advection  of 
momentum  neglected. 

The  observations  of  Lemone  and  Pennell  (1974)  imply 
that  the  vertical  shear  of  the  geostrophic  wind  in  the  lower 
atmosphere  of  the  trade  wind  region  can  be  too  large  to  allow 
vertical  uniformity  of  momentum  in  the  subcloud  mixed  layer. 
Consequently,  we  will  not  assume  that  momentum  is  independent 
of  pressure  in  the  subcloud  mixed  layer.  This  prevents 
integration  of  (29)  and  (30)  over  mixed  layer  depth  and 
requires  computation  of  the  pressure  dependence  of  u and  v. 

Thus,  the  pressure  dependence  of  w'u*  and  w'V' 
must  be  parameterized  in  terms  of  mean  field  quantities. 

To  accomplish  this  we  adopt  an  eddy  diffusion  coefficient 
representation  by  assuming 
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p-g-  K (p)  g 


2 2 * . 3v 

p K K (p)  3^ 


for  p < p. 


where  P1  is  the  pressure  at  the  first  model  gridpoint  above 
the  sea  surface.  Below  level  p^  the  eddy  momentum  flux  is 
determined  from  the  surface  boundary  layer  theory  discussed 
in  Section  2.4. 

We  calculate  the  eddy  diffusion  coefficient, 
K*(p),  from  the  model  of  O’Brien  (1970).  In  this  formula- 
tion the  values  of  the  diffusion  coefficient  and  its  deriva 
tive  in  the  surface  layer  and  at  the  top  of  the  mixed  layer 

are  used  to  derive  a Hermite-interpolating  polynomial  for 

* 

K (p)  given  by 


K*(p)  = K* 


<Pb>  + (^^l)  K*(pi>  - K,(V 


* <p  - Pa>  |p  K*<P)I 


2(K*(Pl) 


- K*<Ph>  /(Ph  - P,> 


where  3K*pP^|pb  is  assumed  zero  and  K*(pb)  is  set  to  a 
small  value  characteristic  of  the  free  atmosphere.  From 
surface  boundary  layer  theory,  K*(pj)  and  — §p—lpi  Rre 
given  by 


r 


The  functional  form  of  the  nondimensional  wind  shear  function, 
4)  £ , is  given  in  Section  2.4.  This  treatment  of  the  eddy 

diffusion  coefficient  has  been  used  by  Pielke  (1973)  and  Lee 
(1974). 
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2.3  THE  CLOUD  LAYER 
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2.3.1  Basic  Equations 


We  now  consider  the  heat  and  moisture  conservation 
equations  (1)  and  (2)  applied  in  the  cloud  layer.  We  assume 
that  w's' , w'q'’ , c,  and  e above  p^  are  due  solely  to  the 
presence  of  cumulus  clouds  that  begin  as  updrafts  originating 
from  the  mixed  layer. 

We  neglect  the  accumulative  storage  of  heat,  water 
vapor,  and  liquid  water  in  the  cloud  ensemble  which  implies 
that  it  is  in  a statistical  steady  state.  This  is  an 
excellent  assumption  in  the  undisturbed  trade  wind  regime 
and  allows  us  to  write  the  balance  of  liquid  water  between 
the  environment  and  the  cloud  ensemble  as 

c - e - |-(^P)o  - R*  = 0 , (33) 


where  R*  is  the  rate  of  raindrop  generation  per  unit  mass 
and  the  subscript  o indicates  that  the  fluxes  are  due  to 
clouds.  Equation  (33)  can  be  used  to  transform  (1)  and  (2) 
into 


If  ♦ 5*5  - = If 
If  * + e || 


Qo  + 


R 7p 


Lv<^>0  - <^>o 


>o  + <^>o 


- R* 


+ Lv**  , 
(34) 

. (35) 
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The  quantities  (uTrFr)o,  (<*>'q')a,  (a and  R*  will  be 
parameterized  in  terms  of  mean  quantities  as  described  in 
Section  2.3.2. 

In  the  cloud  layer,  we  again  assume  that  the  clouds 
do  not  produce  a vertical  flux  of  momentum  and  that  the 
Rossby  number  is  <<1.  Hence,  the  equations  for  conservation 
of  momentum  in  the  cloud  layer  take  the  same  form  as  (29) 


2.3.2  Cumulus  Parameterization 


The  cumulus  parameterization  theory  used  in  this 
model  is  due  to  Arakawa  and  Schubert  (1974).  This  theory 
provides  a closed  parameterization  for  the  interaction  of  a 
cumulus  cloud  ensemble  with  the  large-scale  environment  and 
will  be  reviewed  here. 


2. 3.2.1  Spectral  Representation  of  the  Cumulus  Ensemble 


The  cloud  ensemble  is  spectrally  divided  into 
subensembles  with  all  the  clouds  in  a particular  subensemble 
having  the  same  constant  fractional  entrainment  rate,  X, 
defined  by 


X 


_=£ 

H w0(X,p) 


3 u>c 

Jp 


L(X,p) 


(36) 


where  «a(X,p)  is  the  average  vertical  p-velocity  inside  the 
clouds  in  subensemble  X at  level  p multiplied  times  the 
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fractional  part  of  the  model  domain  covered  by  subensemble 
, and  H = RT/g.  It  is  convenient  to  express  u>0(A,p)  as 

w0(A,p)  - - g n(A,p)  m^A)  dA  , (37) 

where  m^A)  dA  is  the  summation  of  the  instantaneous  cloud- 
induced  vertical  mass  fluxes  at  pb  over  all  the  clouds  in 
subensemble  A and  n(A,p)  is  a nondimensional  normalizing 
function,  equal  to  1 at  pfe.  Equations  (36)  and  (37)  can  be 
combined  to  yield 


! 


( 


■|~(  A ,p)  " - n(A,p)  , (38) 

which  is  a differential  equation  expressing  conservation  of 
mass  for  the  subensemble  that  can  be  solved  for  n(A,p). 

Cloud  tops  are  taken  to  be  the  level  of  vanishing 
buoyancy,  pd(A),  and  the  subensemble  is  assumed  tv  have  a 
thin  layer  of  detrainment  at  this  level  in  which  the 
vertical  mass  flux  goes  to  zero.  Thus,  n(A,p)  = 0 for 
p _<  Pd(A). 


Note  that  since  the  members  of  the  subensemble 
are  assumed  to  be  at  random  phases  in  their  life  cycle, 
mb(A)  dA  is  proportional  to  the  mass  flux  of  a single  cloud 
in  subensemble  A averaged  over  its  entire  lifetime,  with 
the  constant  of  proportionality  equal  to  the  number  of 
clouds  in  the  subensemble.  Note  also  that  because  of  the 
effect  of  entrainment  on  cloud  buoyancy  and  from  the  entrain- 
ment-radius relationship  (see  Simpson,  1971),  small  A 
corresponds  to  deep  clouds  with  large  radii  and  large  A 
corresponds  to  shallow  clouds  with  small  radii  (see  Figure  1). 
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Now,  again  invoking  the  assumption  that  the 
fractional  area  covered  by  active  clouds  is  small,  the  total 
cloud-induced  vertical  fluxes  of  mass,  dry  static  energy, 
water  vapor,  and  liquid  water  averaged  horizontally  across 
the  model  domain  can  be  written  (Schubert,  1974) 


d(p) 

w0(p)  = - g J n(X,p)  mb(X)  dX 


(39) 


Aj(P>  . - i 

(u>'s')0  * - g J n(  X ,p)  |s0(  X ,p)  - s(p)|mb(X)  dX  , (40) 


X (p)  _ . 

(I7q’)0  = - g f n(X,p)  |q0(X,p)  - q(p)|mb(X)  dX  , (41) 


(<^**), 


■ gl 


xd(p) 


n(  X ,p  ) l(X,p)  n^(X)  dX 


(42) 


where  s (X,p).  q (X.p).  and  t(X,p)  are  the  dry  static  energy. 
0 0 

water  vapor  mixing  ratio,  and  liquid  water  mixing  ratio  in- 
side type  X clouds  at  level  p and  Xd(p)  is  the  entrainment 
rate  of  the  clouds  whose  tops  are  at  level  p.  Thus,  the 
total  cloud-induced  flux  of  a particular  quantity  at  level  p 
is  an  integral  of  the  subensemble  fluxes  of  this  quantity 
over  all  the  subensembles  which  penetrate  level  p. 


In  addition,  we  can  express  the  ensemble  average 
raindrop  generation  rate  by 
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(43) 


RS  - / 

o 


xd(P) 


n(A,p)  r(A,p)  m^A)  dX 


where  r(X,p)  is  the  raindrop  generation  rate  for  subensemble 
X at  level  p per  unit  cloud  mass  flux.  We  will  parameterize 
r(A,p)  by 


r(A,p)  = CQ  £(A,p) 


(44) 


1« 


where  CQ  is  a constant.  Furthermore,  we  neglect  the  re- 
evaporation of  raindrops  and  assume  that  all  the  detrained 
liquid  water  evaporates  into  the  environment  immediately. 

To  close  the  cumulus  parameterization  theory, 
sa(A,p),  q0(A,p),  £.(  A , p) , and  m^A)  must  be  determined  as 
a function  of  the  large-scale  quantities.  This  can  be 
achieved  from  the  considerations  which  follow. 


2. 3.2.2  Prediction  of  the  Thermodynamic  Properties 
of  the  Clouds 


The  equations  for  conservation  of  heat,  water 
vapor,  and  liquid  water  in  subensemble  X can  be  written 
(Schubert,  1974) 


fp  bx'*»  so(X*p)l  * 


As(p)  + Ly  c( X , P)|  " 0 


(45) 
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f-  | n(  A ,p)  q0(A,p)|  + % jxq(p)  - c(A,p)|  - 0 

P 

(46) 

] 

. 

_ 

|p  |n(X,p)  UA,p)|  + | c ( X , p ) - r(A,p)|  - 0 ' 

(47) 

The  first  terms  on  the  left-hand  sides  of  (45)  through  (47) 
represent  convergence  of  heat,  moisture,  and  liquid  water, 
respectively,  in  the  subensemble  per  unit  cloud  mass  flux. 

The  terms  in  (45)  and  (46)  which  have  X as  a factor 
represent  entrainment  of  environmental  heat  and  moisture, 
respectively,  per  unit  cloud  mass  flux.  Finally,  the  second 
and  third  terms  in  (47)  represent  the  liquid  water  source 
due  to  condensation  and  the  liquid  water  sink  due  to  raindrop 
formation  in  the  subensemble  per  unit  cloud  mass  flux.  Note 
that  radiational  heating  inside  the  clouds  has  been  neglected, 
which  is  reasonable  since  it  is  small  compared  to  the  other 
terms  in  (45).  Equations  (45)  through  (47)  are  valid  in  the 
region  pb  > p > Pd( A) . 

Equations  (45)  and  (46)  can  be  combined  to  yield 

|n(  A ,p)  h0(A,p)|  + h h(p)  - 0 . (48 ) 

where  h0(A,p)  is  the  moist  static  energy  of  the  clouds 
given  by 

h0(A,p)  ■ s ( A ,p ) ♦ Ly  qc(A,p)  , (49) 
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and  h(p)  is  the  moist  static  energy  of  the  environment  given 


h(p)  = s(p)  + Ly  q(p) 


(50) 


Similarly,  (46)  and  (47)  can  be  combined  to  yield 


n(X,p)  ( q0( A ,p)  + £( A ,p) ] 


H q(p)  - ” plA.PJ  r(  X ,p)  *=  0 
P P 


Furthermore,  above  cloud  base  q0(X,p)  can  be  determined  from 
the  saturation  condition 


q0(A,p)  = q*(p)  + 3-^—  ±-  h0(X,p)  - h*(p)|  , (52) 


where  q*(p)  is  the  saturation  mixing  ratio  of  the  environment 
h*(p)  is  the  saturated  moist  static  energy  of  the  environment 


h*(p)  = s ( p ) + Ly  t»*(p) 


(53) 


y = _x  /ia i) 

Y ‘ Op  'p 


(54)* 


♦(See  Arakawa  and  Schubert,  1974) 
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The  analytical  forms  for  q*  and  X used  in  this  model  are 
given  in  Appendix  B. 

From  our  assumption  that  the  clouds  have  ’’roots’' 
extending  down  into  the  mixed  layer  which  is  consistent 
with  the  observations  of  Pennel  and  Lemone  (1976),  we  take 

so(X’pb>  * sm  * (55) 

q0(A,pb)  « qm  . (56) 

le  further  assume  pc  to  be  at  the  lifting  condensation 
level  and  note  that 

t(X,p)  * 0 for  p > pc  . (57) 


Equation  (44)  along  with  the  differential  equations  (48) 
and  (51)  and  the  boundary  conditions  (55)  through  (57), 
allow  prediction  of  the  thermodynamic  properties  of  the 
cloud  ensemble  in  the  following  manner.  For  pc  < p < pb> 
(48)  and  (51)  are  integrated  upwards  from  pb  taking  into 
account  (55)  through  (57).  From  the  values  of  h0(X,p)  and 
qo(X,p)  so  obtained.  so(X,p)  can  be  determined  from  (49). 
For  pc  > p > Pd(X),  (48)  and  (51)  are  Integrated  upwards 
from  p with  consideration  of  (44).  This,  along  with  (52) 

C 

and  (38)  yields  so(X,p),  qo(X,p),  and  l(X,p).  Cloud  tops 
are  assumed  to  occur  where  the  clouds  lose  buoyancy  which 
Implies 


8vo(Xd(p)’p)  " SV(P)  ’ (58) 
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at  p - Pd(X),  where  sy(p)  is  given  by  (21)  and 

svo(A,p)  - so(A,p)  + cp  t(o .61q  (A,p)  - *(A,p))  . 

(59) 


2.3.2. 3 Prediction  of  m^C  A) 


The  cumulus  parameterization  will  be  closed  if  we 
can  predict  the  cloud-induced  vertical  mass  flux  at  p^, 
mp(A).  To  do  this  we  introduce  the  "cloud  work  function", 
A( A) , given  by 


A(  A) 


Pd(A) 

— / n(A,p)  (svo(A,p) 
P r> 


(60) 


which  is  simply  an  integral  measure  of  the  buoyancy  force 
for  cloud  type  A weighted  with  the  function  n(A,p).  The 
time  derivative  of  A(A)  can  be  written  in  the  form 


dA(  A) 
dt 


/max 

K( A , A ' ) mb(A')  dA 


F(A) 


(61) 


where  K( A , A ' ) and  F( A ) can  both  be  expressed  in  terms  of 
large-scale  quantities  alone. 

The  term  F( A ) in  (61)  is  the  rate  at  which  the 
large-scale  forcing  increases  the  cloud  work  function  for 
subensemble  . The  dominant  effects  in  F( A ) are  typically 
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cooling  above  cloud  base  due  to  infrared  radiation  and 
adiabatic  expansion.  The  kernel,  K(A,A')  is  the  rate  at 
which  type  A'  clouds  increase  the  cloud  work  function  for 
type  A clouds  per  unit  cloud  mass  flux.  It  is  generally 
negative  and  the  dominant  effect  it  represents  is  warming 
above  cloud  base  due  to  adiabatic  compression  associated 
with  compensating  subsidence  in  the  environment  between  the 
clouds.  Hence,  the  first  term  on  the  right-hand  side  of 
(61)  represents  the  rate  at  which  the  cumulus  ensemble 
decreases  the  cloud  work  function  for  subensemble  A through 
release  of  the  moist  convective  instability. 

Arakawa  and  Schubert  (1974)  demonstrate  that  in 
the  tropics  the  time  scale  for  changes  in  the  large-scale 
forcing  is  much  longer  than  the  time  scale  for  the  cloud 
ensemble  to  adjust  to  the  forcing.  This  implies  that  the 
cumulus  ensemble  follows  a sequence  of  quasi-equilibria 
with  the  large-scale  forcing.  Furthermore,  the  cloud  work 
function  in  this  set  of  quasi-equilibria  remains  approxi- 
mately constant  in  time  at  some  relatively  small  value. 

Thus,  (61)  can  be  written 

/max 

J K(A.A')  m^A')  dA'  + F(A)  « 0 . (62) 


The  kernel,  K( A , A ' ) can  be  expressed  as  the  sum 
of  three  parts 


K(  A , A')  « Ky( A , A ' ) + Kd( A , A ' ) + KM( A ) . (63) 
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The  "vertical  mass  flux  kernel",  Ky(X,X'),  accounts  for  the 
decrease  of  the  cloud  work  function  for  type  X clouds  due 
to  compensating  subsidence  in  the  environment  and  is  the 
dominant  effect  in  K(X,X')  as  mentioned  earlier.  The 
"detrainment  kernel",  KD(X,X'),  accounts  for  the  increase 
of  the  cloud  work  function  for  type  X clouds  due  to  cooling 
and  moistening  of  the  environment  caused  by  detrainment  of 
liquid  water  and  saturated  air  from  shallower  clouds.  The 
"mixed  layer  kernel",  Kjj(X)  accounts  for  the  decrease  of  the 
cloud  work  function  for  type  X clouds  due  to  changes  in  the 
thermodynamic  characteristics  and  thickness  of  the  mixed 
layer . 


Analytical  expressions  for  Ky(X.X'),  KD(X,X'), 
Km(X),  and  F(X)  can  be  obtained  by  formally  taking  the  time 
derivative  of  (60)  and  making  use  of  earlier  expressions. 

The  details  of  this  calculation  are  given  in  Arakawa  and 
Schubert  (1974)  and,  because  of  their  great  complexity, 
will  not  be  presented  here. 

For  use  of  the  theory  in  a numerical  prediction 
model,  it  is  unnecessary  to  evaluate  an  analytical  expression 
for  F(X).  Rather,  it  is  convenient  to  replace  F(X)  by  the 
equivalent  expression 


F(  X ) 


ALS^ 

At 


where  ALg(X)  is  the  cloud  work  function  produced  by  the 
large-scale  forcing  alone  during  one  time  step,  At,  of 
the  model.  Hence,  (62)  can  be  written 


( 
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X v Aw  „(  X) 

G(  X)  - /max  K(X.X')  m^X)  dX'  + -L~  - 
Jo 

(64) 

Using  (60),  (63),  and  the  expressions  in  Appendix  B of 
Ar&kawa  and  Schubert  (1974),  equation  (64)  (which  is  a 
Fredholm  integral  equation  of  the  first  kind)  is  solved  for 
mb(X),  subject  to  the  constraints 


mb(X)  = 0 

for 

A(  X ) < 0 

or 

o 

V 

r< 

o 

mb(X)  > 0 

for 

G( X ) > 0 

. 

(65) 

Additional  changes  in  s and  q are  then  calculated  from  the 
cloud  terms  in  (36)  and  (37).  The  integral  equation,  (64), 
is  solved  with  a simple  iterative  scheme  described  in 
Appendix  C which  is  similar  to  that  discussed  by  Schubert 
(1973). 


Note  that  this  cumulus  parameterization  theory 
does  not  explicitly  predict  o.  Following  Ogura  and  Cho 
(1973),  however,  we  can  estimate  o by  assuming  that  the 
magnitude  of  the  vertical  velocity  of  the  updrafts  entering 
the  clouds  at  pb,  w^,  is  independent  of  X and  o.  Thus  we 
have 


Pb  B wb 


(66) 
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As  an  estimate  of  the  fractional  cloud  coverage,  where  pfe 
f is  the  density  of  the  air  at  and  is  calculated 

from  (39).  Ogura  and  Cho  found  wfa  to  be  about  1 m s”1 
from  their  diagnostic  calculations,  which  is  the  value  used 
in  this  model, 
f 


2. 3.2.4  Detrainment  Form  of  the  Cumulus  Parameterization 
Theory 

It  is  instructive  to  write  the  cloud  terms  of 
(34)  and  (35)  in  "detrainment  form"  (see  Schubert,  1974). 
Differentiating  (40)  through  (42)  with  respect  to  p and 
making  use  of  (38)  and  (45)  through  (47),  the  conservation 
equations  for  heat  and  moisture,  (34)  and  (35),  can  be 
written 

g ♦ r*i  ♦ =g  - ♦ D<P>I=  - 5 - “I  + “o(p)  If  • 

(67) 

+ ^*Vq  + = D(p)|q  - q + t|  + Q( P ) , (68) 


where  the  D(p) 
given  by 

D(p) 

and 


The  total  detrainment  rate  is  ^ero  below  the  lowest  detrain 
ment  level.  Thus,  below  this  level  the  ensemble  modifies 
the  large-scale  fields  simply  through  cloud-induced  subsi- 
dence in  the  environment  as  represented  by  the  terms 
w_(p)  U and  w (p)  |a  in  (67)  and  (68). 
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2.4  THE  SURFACE  LAYER 


Under  steady,  convectively  unstable,  horizontally 
homogeneous  conditions,  Monin-Obukhov  similarity  theory 
predicts  that  a surface  layer  of  thickness  -L  exists  in 
which  the  eddy  fluxes  depart  little  from  their  surface 
values  and  properly  scaled  derivatives  of  temperature, 
moisture  and  momentum  are  universal  functions  of  z/L . 

Here,  z is  the  height  above  the  surface  and  L is  the  Monin- 
Obukhov  length  defined  by  (see  Kraus,  1972) 


L 


pa  u* 


c T-  + O.eiC^Tq  )o 
P O 


(71) 


where  u*  is  the  surface  friction  velocity  given  by 


(72) 


( 


€ 


and  is  the  surface  wind  stress  vector.  Other  symbols 
o 

are  defined  in  Appendix  A.  Note  that  L is  negative  in  a 
convectively  unstable  boundary  layer  which  is  the  case  we 
restrict  our  attention  to  here. 

Thus,  in  the  surface  layer  we  have 


32  " Oa  8 KU*  2 Cp 


(73) 


1 
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ll  - (u>  q >o  *q(z) 

3Z  Pa  g KU*  Z 


(74) 


3U  _ u*  A /z\ 

75  « V(l) 


(75) 


where  U is  wind  speed.  Ths  forms  of  the  4>  functions  have 
been  reviewed  by  Hogstrom  (1974)  and  it  is  generally 
accepted  that  most  atmospheric  data  are  well  represented  by 


which  are  the  forms  used  here. 


Using  (76).  the  profile  relations  (73)  through 
(75)  can  be  integrated  from  zq , the  surface  roughness 
length,  to  z (see  Paulson,  1970) 


(77) 
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"(I)  ■ qo  * P,  g~u.  ln| 


2 In 


(^)l 


(78) 


U(z)  = UQ  + ln(f-)  - 2 ln(— § — ) 
ln( S^1-)  + 2 tan_1  ym  " iTr  ’ 


where  6o>  qQ,  and  UQ  are  potential  temperature,  water  vapor 
mixing  ratio,  and  wind  speed  at  zQ , and 


*e  ■ >'q 


(>  - 9 r) 


- (!  - 15  e) 


(80) 


(81) 


At  a height  of  -L,  6 and  q are  given,  to  excellent  approxi- 
mation, by  their  mixed  layer  values.  Thus,  setting  z * -L 
in  (77)  and  (78)  yields,  after  rearrangement 

_ (em  * eo) 

<“  * >o  ■ S c.  8 TT773\ 1 

lnU_  - 1.466 


- B H, 
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(82) 


and 


- *cu,  p g 


(%  - %) 


In 


m- 


1.466 


g Q, 


(83) 


where  H is  the  surface  sensible  heat  flux  and  L is  the 

O VO 


surface  latent  heat  flux  (both  taken  positive  downward). 
Since  we  do  not  assume  that  momentum  is  uniform  in  the  mixed 
layer  it  is  not  advantageous  to  integrate  (75)  to  the  top  of 
the  surface  layer.  Rather,  we  integrate  (75)  from  zQ  to  z1 
where  is  the  height  of  the  first  model  gridpoint  and 


z^  < -L. 


Thus,  (79)  yields  after  rearrangement 


(° 


:1  - «o)/ [ 1 n i;  - 2 '"((-ir1)  - lB( 


♦ 2 tan_1y1  - in  , (84) 


where  is  the  wind  speed  at  the  first  model  gridpoint, 
UQ  is  the  wind  speed  at  zQ,  and 


Following  Clarke  (1970)  and  others,  we  assume 
that  the  roughness  length  for  the  sea  surface  is  given  by 

zQ  * max  (o. 032  u*2/g,  zmin)  . (86) 

where  zmin  * 0.0015  cm.  We  set  6q  equal  .to  the  sea  surface 
temperature  and  q equal  to  0.98  times  the  saturation  value 
at  the  sea  surface  to  account  for  the  effect  of  salinity 
on  the  vapor  pressure  of  water  (see  Roll,  1971).  We  set 
U equal  to  the  surface  drift  velocity  which  is  approximated 
by  (Kraus, 1977) 


i'o 


v + 
-op 


Xc 


(87) 


where  p and  p are  reference  densities  for  air  and  water, 

Si  ^ 

v is  the  mean  current  vector  in  the  upper  layer  of  the 
-*-og 

ocean.  The  direction  of  the  surface  wind  stress  vector  is 
assumed  to  be  that  of  the  wind  at  . The  wind  vector  at 
is  predicted  by  (29)  and  (30).  Thus,  using  the  fact  that 


which  follows  from  (6),  (12),  and  (13),  equations  (82) 
through  (88)  allow  prediction  of  the  surface  fluxes  of  heat, 
moisture,  and  momentum  in  terms  of  model -predicted  quantities. 
We  solve  (71),  (72)  and  (82)  through  (86)  iteratively  at  each 
time  step  and  assume 
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2.5  RADIATIVE  TRANSFER  EQUATIONS 


The  net  radiational  heating  rate  of  the  atmosphere 
averaged  across  the  model  domain  can  be  written 


QR  - - g 


(1  " o)  3p(F£lR  + S£lr) 


° 3p(F£LD  + S£ld)| 


(91) 


where  F£LR  and  S*LR  are  the  net  downward  fluxes  of  long- 
and  short-wave  radiation  in  the  clear  region,  and  F^^  and 
S£ld  are  the  net  downward  fluxes  in  the  cloudy  region.  Note 
that  since  a <<  1 for  the  physical  system  to  be  modeled,  the 
net  radiational  heating  due  to  the  presence  of  clouds  is  very 
small.  Thus,  (91)  can  be  written  to  excellent  approximation 
as 


QR  = " g 3p(F£lR  + S£lr)  . (92) 


We  calculate  the  long-  and  short-wave  fluxes  in  the  manner 
described  in  the  next  two  sections. 


2.5.1  Long-Wave  Radiative  Flux 

Using  the  approximations  of  Danard  (1969),  we 
calculate  the  terrestrial  long-wave  radiative  flux  at  level 
p in  the  clear  region  from 
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f*(p)clr  * 


- e Bfl  + 
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| [£  f4<P'»rtP.P*)  5(p'>  dp 

rP o / \ / ^\0 • 85 / T \0*5 

^ (c0To  - rl(p',)E*<p’p')  q(p')  (§;)  dp 


for  p < p < p 


(93) 


where  B is  the  Stefan-Boltzmann  constant,  eo  is  the  emis- 
sivity  of  the  sea  surface,  Tr  is  a reference  temperature, 
Pr  is  a reference  pressure,  pe  is  the  pressure  at  the  top 
of  the  model , and 


e*(p,p')  * e'(|w(p)  - w(p')|) 


Here,  w(p)  is  the  pressure  corrected  precipitable  water 
vapor  defined  by 


(94) 


and 
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We  assume  that  the  precipitable  water  at  pe  is  a constant 
and  we  approximate  e'  by  fitting  Elsasser's  (1960)  experi- 
mental curve  for  the  flux  emissivity  of  pure  water  vapor, 
e,  with  an  analytical  function. 

The  first  term  on  the  right-hand  side  of  (93)  is 
the  black  body  emission  from  the  sea  surface.  The  first 
and  second  integrals  on  the  right-hand  side  of  (93) 
represent  the  contribution  to  the  long-wave  flux  at  level 
p due  to  emission  and  absorption  of  radiation  above  and 
below  level  p respectively.  See  Haltiner  (1971)  for  a 
detailed  discussion  of  this  radiative  transfer  model. 

Note  that,  although  (92)  and  (93)  give  a good 
prediction  of  QR,  this  long-wave  radiative  transfer  model 
overestimates  the  infrared  cooling  of  the  sea  surface. 
Consequently  we  use  the  empirical  relationship  of  Reed 
(1976)  to  calculate  the  net  long-wave  flux  at  the  sea 
surface.  Thus 

F*  = T4  0.254  - 0.811  q,ft  , (96) 

O O O 10 

where  q^Q  is  the  water  vapor  mixing  ratio  at  10  m height 
and  is  calculated  from  the  profile  relation  (78). 

2.5.2  Short-Wave  Radiative  Flux 


Following  Arakawa,  Katayama,  and  Mintz  (1968), 
we  assume  the  short-wave  radiation  to  consist  of  two  parts: 
(a)  with  wavelength  less  than  0.9  y in  which  Rayleigh 
scattering  is  significant  but  absorption  by  water  vapor  is 
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negligible,  and  (b)  with  wavelength  greater  than  0.9  u in 
which  absorption  by  water  vapor  is  significant  but 
scattering  is  slight.  In  this  formulation  the  incident 
radiation  is  partitioned  according  to 

SS  * 0.651  S cos  a 
o o 

S*  - 0.349  So  cos  a , (97) 


where  S is  the  net  incident  solar  radiation  at  the  model 

e® 

top,  S is  the  downward  flux  of  solar  radiation  subject  to 
° A 

scattering,  S"  is  the  downward  flux  of  solar  radiation 
subject  to  absorption,  and  a is  the  Zenith  angle.  Then, 
Si.p(p)  is  represented  by  the  empirical  relationship 


SWP>  - 


+ Sk 


1 - C (w*(p)  sec  a) 


0.303 


1 - ao 
1 - aoas 


for  p < p < p„ 
*e  — — o 


(98) 


where  a is  the  albedo  of  the  scattered  radiation  and  is 
o 

taken  to  be 


0.085  - 0.245  log 


(Lo  cos  a) 


and 

w*(p)  - w(pc)  " W(P> 


(99) 


(100) 
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C is  a constant,  and  a„  is  the  albedo  of  the  sea  surface, 
s s 

The  expression  (1  - a a ) in  the  denominator  of  the  second 

o s 

term  on  the  right-hand  side  of  (98)  is  present  to  model  the 
enhancement  of  the  downward  flux  of  radiation  due  to  back 
scattering  from  multiple  reflections.  This  short-wave 
radiative  transfer  model  is  described  in  more  detail  by 
Haltiner  (1971). 

In  order  to  determine  the  daily-averaged  value 
of  the  solar  heating  rather  than  the  instantaneous  value, 

SQ  cos  a is  averaged  over  24  hours  with  the  constraint 
S cos  a = 0 for  a > $u  and  a < -in.  This  yields 


/S  cos  o')  = 

\ o / o 


(sin  <fr  sin  6) 


, (cos  $ cos  6) 

2tt 


f cos  ( d ( 


(101) 


where 

4>  * latitude  at  center  of  model  domain 
6 * declination  angle  of  sun 
0 * cos-1  (tan  4 tan  d)  , 

and  < > indicates  the  average  over  24  hours.  We  use 
<SQ  cos  a>  in  (97)  through  (99)  and,  although  this  is  not 
exactly  equivalent  to  averaging  (97)  through  (99)  over  24 
hours,  it  is  adequate  for  our  purposes. 
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Finally,  the  net  short-wave  flux  at  the  sea 
surface  is  calculated  from 

SS  - SCLR(Po) 

where  again  we  have  made  use  of  the  fact  that  o <<  1 for 
the  cases  to  be  considered  here. 


(102) 


2.6  ONE-DIMENSIONAL  REPRESENTATION 


We  will  adopt  a one-dimensional  model  by  inte- 
grating equations  (26)  through  (30)  and  (34)  and  (35)  with 
^•Vq,  ^b*^Pb»  ug,  vg,  «.  pQ,  and  the  sea  surface  temperature 
imposed.  Note  that  specification  of  the  vertical  profile  of 
the  geostrophic  wind  implies  Vs  in  the  following  manner. 

f\  - - *(!f)p  • <103) 

' - «(H)p  • (104) 

and  the  horizontally  averaged  thermal  wind  equations  with 
* the  effect  of  moisture  on  density  neglected 


(105) 


(106) 


the  horizontally  averaged  gradient  of  dry  static  energy  on 
a level  of  constant  pressure  can  be  written 
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(107) 


where  i and  j are  unit  vectors  in  the  x and  y directions. 
Hence,  the  horizontally  averaged  horizontal  advection  of 
dry  static  energy  can  be  written 


+ 


+ 


(108) 


Thus,  the  horizontal  advection  of  heat  is  included  in  a self- 
consistent  fashion. 

For  all  of  the  experiments  reported  here  we  will 
set  v*7q  equal  to  zero.  This  is  a fairly  good  assumption 
for  the  undisturbed  trade  wind  regime.  For  example, 

Esbensen  (1975)  analyzed  the  BOMEX  data  and  found  that  the 
eddy  flux  of  moisture  at  the  top  of  the  mixed  layer  nearly 
balanced  the  sea  surface  evaporation  during  the  period  of 
his  study.  This  along  with  the  well-known  steadiness  of 
the  trade  wind  thermodynami c structure,  implies  that  the 
horizontal  advection  of  moisture  in  the  mixed  layer  is,  on 
the  average,  relatively  small  compared  to  the  other 
processes  that  contribute  to  the  mixed  layer  moisture 
balance . 
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In  addition,  we  set  ec*ua*  *°  zero  in  ftii 

of  the  experiments  reported  here.  This  is  a good  assumption 
since  the  mixed  layer  thickness  tends  to  be  spatially  uni- 
form in  the  undisturbed  trade  wind  regime.  Consequently, 
the  contribution  of  ^b*^Pb  in  (28)  is  expected  to  be 
negligible  compared  to  the  other  terms. 
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2.7  COMPUTATIONAL  CONSIDERATIONS 


2.7.1  Model  Grid 


In  order  to  resolve  the  surface  layer  without 
demanding  unnecessary  core  storage,  a vertically  stretched 
grid  is  employed  with  the  stretching  function  given  by 


1 - b-  j2  - b, 

l ‘ 


with 


(109) 


bj  = 1.106  x 10 


and 

b2  - 1.175  x 10 


Values  of  s,  q,  u,  and  v are  computed  at  levels  defined  by 
setting  j * 1,2, 3,..., 28  in  (109).  This  choice  of  para- 
meters places  the  first  gridpoint  at  approximately  20  m 
above  the  sea  surface  and  the  last  gridpoint  at  about 
tropopaule  level.  The  upper  boundary  conditions  on  the 
prognostic  variables  are  given  in  Table  1. 

Note  that  pb  as  determined  by  (28)  does  not  fall 
at  a model  gridpoint  in  general.  Hence,  in  order  to 
calculate  As  and  Aq , we  utilize  a "floating  gridpoint"  at 
Pg  which  moves  up  and  down  with  the  top  of  the  mixed  layer. 
Assuming  no  clouds  detrain  at  pg,  we  can  write 
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TABLE  1.  (Continued) 
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TABLE  1.  (Conti 
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With  and  q (p^ ) predicted  by  (110)  and  (111),  As  and 

Aq  are  calculated  from  (16). 


2.7.2 


Numerical  Scheme 


We  adopt  a five  cloud-type  cumulus  ensemble  by 
choosing  the  values  of  X^  listed  in  Table  1.  All  integrals 
over  X as  well  as  all  integrals  over  pressure,  are  per- 
formed by  using  the  trapazoidal  rule.  The  differential 
equations  (48)  and  (51)  are  solved  with  a modified  Euler 
predictor-corrector  scheme. 
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Below  the  lowest  level  of  detraimnent , (67)  and 
(68)  are  used  in  place  of  the  equivalent  expressions  (34) 
and  (35).  This  is  done  because  the  detrainment  form  of 
the  equations,  (67)  and  (68),  appear  to  have  computational 
advantages  over  (35)  and  (36)  at  the  top  of  the  transition 
layer  and  near  cloud  base.  The  finite  difference  scheme  is 
forward  in  time  with  the  eddy  diffusion  terms  in  (29)  and 
(30)  treated  implicitly  and  all  the  vertical  advective  terms 
in  (29),  (30),  (34),  (35),  (67),  (68),  (110),  and  (111) 
computed  with  upstream  differences. 


2.7.3  Large-Scale  Condensation  and  Dry  Convective 
Adjustment 


The  cumulus  parameterization  theory  described  in 
Section  2.3.2  accounts  for  subgrid  scale  condensation, 
evaporation,  and  precipitation  which  can  occur  when  the 
atmosphere  is  not  saturated  in  the  large-scale  sense.  For 
generality,  however,  we  must  further  account  for  large- 
scale  moisture  phase  changes  when  the  atmosphere  becomes 
saturated  across  the  entire  model  domain.  To  do  this,  the 
adjustment  scheme  of  Hack  and  Schubert  (1976)  is  implemented 
in  the  model.  It  turns  out,  however,  that  for  the  experi- 
ments reported  here  large-scale  saturation  is  never  achieved. 

Similarly,  if  the  potential  temperature  decreases 
with  altitude  in  some  region  above  pfc  during  the  course  of 
the  integration,  subgrid  scale  dry  convection  is  assumed 
to  occur  such  that  a dry  adiabatic  lapse  rate  results 
there.  This  is  accomplished  by  use  of  a dry  convective 
adjustment  scheme  described  by  Hack  and  Schubert  (1976). 
However,  for  the  experiments  described  here,  an  unstable 
stratification  requiring  dry  convective  adjustment  never 
occurs. 
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2.8 


SUMMARY  OF  MODEL  FORMULATION 


A one-dimensional  model  of  the  tropical  marine 
boundary  layer  has  been  proposed  and,  for  convenience,  will 
be  summarized  here. 


2.8.1  The  Surface  Layer 


A surface  layer  extending  from  the  sea  surface 

to  a level  indicated  by  p is  assumed  to  exist  in  which  the 

a 

Monin-Obukhov  Similarity  Theory  applies.  This  allows 
determination  of  the  surface  fluxes  in  terms  of  model- 
predicted  quantities. 


2.8.2  The  Subcloud  Mixed  Layer 


A subcloud  mixed  layer  is  assumed  to  exist 
extending  from  pa  to  a level  denoted  by  pb  in  which  s and 
q but  not  u and  v are  independent  of  altitude.  Conservation 
equations  for  heat,  moisture,  mass,  and  momentum  in  this 
region  are  (26)  through  (30). 

The  subcloud  mixed  layer  is  topped  by  a stable 
transition  layer  which  is  modeled  as  a discontinuity  in 
s and  q. 
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2.8.3  The  Cloud  Layer 

A cloud  layer  in  which  s,  q,  u,  and  v are  all 
allowed  to  vary  with  altitude  exists  extending  from  pfa  to 
the  model  top.  All  vertical  subgrid  scale  fluxes  of  heat  i 

and  moisture  and  all  subgrid  scale  condensation  and 
evaporation  in  this  region  are  assumed  to  be  accomplished 
by  cumulus  clouds  that  begin  as  updrafts  originating  from 
the  mixed  layer  and  have  their  bases  just  above  pb>  These 
clouds  are  parameterized  using  the  theory  of  Arakawa  and 
Schubert  (1974).  Conservation  equations  for  heat  and 
moisture  in  this  region  are  (34)  and  (35).  The  equations 
for  the  conservation  of  momentum  in  this  region  take  the 
same  form  as  (29). 

2.8.4  Large-Scale  Forcing  Parameters 

Because  this  is  a one-dimensional  study,  o>,  ug. 
vgI  ^*Vs,  j*Vq,  and  ]£b*Vpb  must  be  imP°sed*  In  addition, 
since  there  is  no  upper  ocean  model  involved,  the  SST  must 
also  be  specified. 

For  all  of  the  experiments  reported  here,  we  set 
^•Vq  * 0 and  j»Vpb  * 0.  The  thermal  wind  equation  predicts 
Vs  as  a function  of  ug  and  vg,  and  thus  the  imposed 
quantities  which  are  regarded  as  large-scale  forcing  para- 
meters in  our  experiments  are  u,  ug,  vg,  and  the  SST. 


2-50 


Section  3 


DESCRIPTION  OF  THE  EXPERIMENTS 


We  will  model  the  undisturbed  trade  wind  regime 
by  finding  quasi-steady  solutions  to  the  equations  for 
various  combinations  of  the  imposed  large-scale  forcing 
parameters.  The  trade  winds  are  particularly  suited  for 
this  type  of  study  because  of  their  well-known  steadiness. 

Three  sets  of  experiments  are  performed.  Each 
set  is  accomplished  by  imposing  w,  ug,  vg,  and  the  SST  and 
integrating  (26)  through  (30),  (34)  and  (35)  from 
arbitrary  initial  conditions  to  a quasi-steady  state  in 
which  the  thermodynamic  variables  change  by  only  about 
0.1%  per  hour.  Then,  with  the  same  values  for  all  the 
imposed  parameters  except  the  SST,  the  integration  is 
repeated  and  other  balanced  states  corresponding  to 
various  SST  perturbations  are  found.  The  resulting 
differences  in  the  thermodynamic  structure  of  the  atmos- 
phere and  the  magnitude  of  the  surface  fluxes  are  then 
examined. 


Five  experiments  are  performed  in  each  set  with 
the  SST  varying  by  0.25°  K in  the  range  from  301.5  to 
302.5°  K.  The  vertical  profiles  of  u>  and  ug  in  the  lowest 
3 km  for  all  three  sets  of  experiments  are  shown  in 
Figure  2 and  Figure  3.  The  N-S  component  of  the  geo- 
strophic  wind  is  assumed  to  be  independent  of  altitude 
and  given  by: 


Experimental  Set 


2 
3 

Model  constants  for  all  experiments  are  given  in  Table  1. 

The  imposed  parameters  for  these  experiments  are 
fairly  representative  of  the  average  situation  in  the  un- 
disturbed summertime  Western  Atlantic  Trades.  Note  also 
that,  consistent  with  the  assumed  profile  of  u^  and  the 
thermal  wind  equation,  increasing  |v  I increases  the 
cooling  rate  in  the  lower  atmosphere  due  to  horizontal 
advection.  Thus,  the  advective  cooling  of  the  marine 
boundary  layer  increases  from  experimental  set  1 to  set  3. 
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Section  4 
RESULTS 

Figure  4 shows  the  vertical  profile  of  potential 
temperature  resulting  from  an  integration  in  experimental 
set  1 with  a SST  of  302.0  °K.  The  layered  structure  of 
the  marine  boundary  layer  discussed  earlier  is  clearly 
evident  with  the  trade  inversion  at  about  1700  m hiehgt 
and  the  top  of  the  mixed  layer  at  about  500  m height. 
Consistent  with  the  diagnostic  studies  alluded  to  in 
Section  1,  the  clouds  are  confined  below  the  top  of  the 
trade  inversion  and  the  c^mpressional  heating  in  the 
inversion  due  to  large-scale  subsidence  is  balanced  by 
the  evaporative  cooling  of  detrained  cloud  droplets.  In 
the  cloud  layer  below  the  base  of  the  inversion,  the  main 
balance  is  between  the  cooling  due  to  radiation  and 
horizontal  advection  and  the  warming  due  to  cumulus- 
induced  subsidence.  In  the  mixed  layer  the  warming 
caused  by  the  surface  and  entrainment  fluxes  of  heat  is 
balanced  by  the  cooling  due  to  radiation  and  horizontal 
advection . 

Figure  5 shows  the  vertical  profile  of  water 
vapor  mixing  ratio  for  the  same  experiment.  The  layered 
structure  is  again  evident  and  in  this  case  the  drying 
effect  of  the  large-scale  subsidence  in  the  trade 
inversion  is  balanced  by  the  detrainment  of  both  liquid 
water  and  water  vapor.  Note  that  Figure  5 indicates 
3q/3p  * 0 in  the  lower  part  of  the  cloud  layer.  This  is 
because  all  of  the  clouds  h»/e  their  tops  in  the  trade 
inversion.  Thus,  the  detrainment  rate,  D(p),  is  zero  in 
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Figure  4. 
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Vertical  profile  of  potential 
temperature  in  the  lower 
atmosphere  for  experimental 
set  1,  SST  = 302.0  K. 
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Figure  5.  Vertical  profile  of  water 
vapor  mixing  ratio  in  the 
lower  atmosphere  for 
experimental  set  1 , 

SST  * 302.0  K. 
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the  lower  part  of  the  cloud  layer,  which  from  (68)  requires 
3q/3p  * 0 there  in  the  steady  state  since  ^*Vq  is  assumed 
zero.  In  the  mixed  layer,  the  surface  evaporation  is 
balanced  by  the  entrainment  flux  at  p^. 

Figure  6 shows  the  vertical  profiles  of  u and  v 
for  the  same  experiment.  In  the  free  atmosphere  the  wind 
remains  near  its  geostrophic  value  while  an  Ekman  balance 
exists  in  the  mixed  layer  where  the  eddy  viscosity  is 
significant.  Thus  in  the  mixed  layer,  as  expected,  the 
wind  vector  is  smaller  in  magnitude  than,  and  directed  to 
the  left  of,  the  geostrophic  wind  vector. 

Consistent  with  the  diagnostic  studies  mentioned 
in  Section  1,  the  cumulus-induced  subsidence  in  the 
environment  at  the  top  of  the  mixed  layer,  ~(w0)b  is  more 
than  an  order  of  magnitude  larger  than  the  large-scale 
subsidence  there,  (see  equation  (8)).  Thus,  in  the 
steady  state,  the  deepening  rate  of  the  mixed  layer  caused 
by  entrainment  at  its  top  is  balanced  essentially  by  this 
cumulus  induced  subsidence. 

The  physical  balances  responsible  for  maintaining 
the  vertical  structures  of  the  marine  boundary  layers  as 
well  as  the  vertical  structures  themselves  for  the  other 
experiments,  are  qualitatively  similar  to  the  case  Just 
described  and  will  not  be  discussed  here.  Instead,  we 
will  examine  relationships  among  the  cloud  mass  flux,  the 
mixed  layer  thermodynamic  variables,  the  surface  fluxes, 
and  the  sea  surface  temperature. 

Figure  7 shows  the  percent  cloud  cover,  assumed 
proportional  to  the  total  cloud  base  mass  flux  according 
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Figure  6. 


Vertical  profiles  of  u( 

and  v( ) in  the  lower 

atmosphere  for  experimental 
set  1,  SST  « 302.0  K. 


to  (66),  as  a function  of  SST  for  the  three  sets  of 
experiments.  As  expected,  increasing  the  SST  increases 
the  cloud  cover.  Furthermore,  the  cloud  cover  is  a weakly 
nonlinear  function  of  SST  with  the  rate  of  increase 
becoming  larger  at  higher  SST's.  Figure  7 also  indicates 
that  the  cloud  cover  at  a particular  SST  increases  from 
set  1 to  set  3.  This  is  because  more  vigorous  cumulus 
convection  is  required  to  balance  the  greater  advective 
cooling  in  the  cloud  layer  as  the  southward  somponent  of 
the  geostrophic  wind  is  made  more  negative. 


Figure  8 shows  mixed  layer  portential  temperature 
for  the  three  sets  of  experiments  as  a function  of  SST. 

Note  that  it  shows  little  variation  from  set  1 to  set  3 
and  increases  linearly  with  the  SST  in  all  cases  with  a 
slope  of  about  1.1.  Figure  9.  on  the  other  hand,  indicates 
a quite  different  behavior  for  the  mixed  layer  water  vapor 
mixing  ratio  which  decreases  significantly  as  |vg|  in- 
creases and  shows  a weakly  nonlinear  decrease  with  SST. 

The  results  indicate,  however,  that  the  average  moisture 
content  of  the  entire  marine  boundary  layer  increases  with 
the  SST,  implying  that  the  cloud  layer  becomes  more  moist 
as  the  SST  increases.  Furthermore,  this  increase  is 
approximately  linear  with  a slope  of  about  0.5  g kg  * °K  *. 


The  net  heat  flux  at  the  sea  surface,  Eq,  can  be 

written 


H + L O + F*  + 
O VO  o 


ss 


where  H is  the  surface  sensible  heat  flux  (calculated 
o 

from  (82)),  LyQ0  the  surface  latent  heat  flux  (calculated 
from  (83)),  F*  the  surface  long-wave  radiative  flux 
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Figure  8. 


SST ( °K ) 


Mixed  layer  potential  temperature  as 
a function  of  SST  for  (a)  experimental 
set  1 ( — - — ).  (b)  experimental  set 

2 ( ),  and  (c)  experimental  set 

3 ( ). 


(calculated  from  (96)),  and  S*  the  surface  short-wave 

radiative  flux  (calculated  from  (102)).  Figures  10 

through  12  show  these  quantities  as  a function  of  SST  for 

the  three  sets  of  experiments.  Note  that  H , F* , and  S* 

o o o 

show  almost  no  variation  from  set  1 to  set  3 and  are 

practically  independent  of  SST  for  the  range  considered. 

The  magnitude  of  the  surface  latent  heat  flux,  however, 

increases  from  set  1 to  set  3 and  varies  linearly  with  the 

SST  for  all  three  sets  with  a slope  of  approximately 

-105  cal  cnT2  day-1  °K~1.  Thus,  essentially  all  of  the 

variation  of  E is  due  to  the  latent  heat  flux  and  we  find 
o 

dE  /dT  about  1.6  times  the  value  estimated  by  Bathan 
O'  o 

et  al . (1976)  and  about  1.1  times  the  value  estimated  by 
Piacsek  et  al.  (1976). 

Figure  13  shows  the  surface  latent  heat  flux 
plotted  against  percent  cloud  cover  for  an  SST  of  302.5°  K. 
The  increase  of  the  surface  evaporation  rate  with  increasing 
cloud  cover  indicated  by  Figure  13  is  consistent  with  the 
findings  of  Ogura  and  Cho  (1973)  who  studies  the  inter- 
action between  the  subcloud  and  cloud  layers  with  a simple 
analytical  model.  This  relationship  between  cloud  activity 
and  surface  evaporation  rate  implies  that  cumulus  convection 
must  be  treated  adequately  in  any  model  used  to  predict  the 
surface  latent  heat  flux  since  feedback  between  the  surface 
evaporation  and  the  cumulus  convection  can  amplify  dQo/dTQ. 
The  neglect  of  cumulus  convection  by  Clancy  (1977)  may  be 
the  reason  for  the  unreasonably  small  value  of  dEQ/dTo 
predicted  by  that  study. 
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Surface  heat  fluxes  (positive  down- 
ward) as  a function  of  SST  for 
experimental  set  1.  The  total  flux 
En  is  the  algebraic  sum  of  the  solar 
flux  Si,  the  long-wave  flux  Fj,  the 
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evaporative  heat  flux. 
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Figure  11.  Surface  fluxes  (positive  downward) 
as  a function  of  SST  for  experi- 
mental set  2. 
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Figure  12.  Surface  fluxes  (positive  downward) 
as  a function  of  SST  for  experi- 
mental set  3. 
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Section  5 

SUMMARY  AND  DISCUSSION 


A one-dimensional  numerical  model  of  the  tropical 
atmosphere  has  been  constructed  and  used  to  investigate  the 
effect  of  sea  surface  temperature  perturbations  in  the 
range  likely  to  be  generated  by  Ocean  Thermal  Power  Plants 
on  the  marine  boundary  layer  of  the  summertime  Western 
Atlantic  trade  wind  regime.  Conservation  equations  for  heat, 
moisture,  and  momentum  are  solved  on  a vertically  stretched 
grid  extending  from  the  sea  surface  to  about  tropopause 
level.  Vertical  eddy  fluxes  of  heat  and  moisture  above  the 
surface  layer  are  accounted  for  by  the  theory  of  Arakawa  and 
Schubert  (1974).  In  the  subcloud  mixed  layer  the  vertical 
eddy  flux  of  momentum  is  modeled  with  the  eddy  diffusion 
coefficient  representation  of  O'Brien  (1970).  In  the 
surface  layer  the  eddy  fluxes  of  heat,  moisture,  and 
momentum  are  determined  from  the  Monin-Obukhov  similarity 
theory.  The  flux  of  long-wave  radiation  is  treated  with 
the  formulation  of  Danard  (1969),  while  the  short-wave 
radiative  flux  is  calculated  from  the  model  of  Arakawa, 
Katayama,  and  Mintz  (1968). 

Vertical  profiles  of  the  large-scale  vertical 
p-velocity  and  the  geostrophic  wind  are  imposed  along  with 
the  SST , surface  pressure,  and  boundary  conditions  for  the 
prognostic  variables.  The  horizontal  gradient  of  moisture 
is  assumed  zero  while  the  horizontal  gradient  of  temperature 
is  implied  from  the  imposed  geostrophic  wind  profile.  Thus, 
the  horizontal  advection  of  moisture  is  ignored  but  the 
horizontal  advection  of  heat  is  included  in  a self- 
consistent  fashion. 
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The  model  is  integrated  from  arbitrary  initial 
conditions  to  a quasi-steady  state.  Then,  for  various 
combinations  of  the  imposed  parameters,  the  integration 
is  repeated  and  other  quasi-steady  states  are  found.  The 
influences  of  SST  and  cloud  cover  (which  is  assumed  pro- 
portional to  the  total  cloud  base  mass  flux)  on  the 
structure  of  the  marine  boundary  layer  and  the  magnitude 
of  the  surface  fluxes  are  then  examined. 

The  height  of  the  trade  wind  inversion  and  the 
thickness  of  the  atmospheric  mixed  layer  are  not  appreciably 
affected  by  changes  of  SST  of  1 °C.  In  addition,  the 
surface  fluxes  of  long-wave  radiation,  short-wave  radiation, 
and  sensible  heat  are  found  to  be  little  affected  by  small 
changes  in  SST  and  cloud  cover.  The  surface  latent  heat 
flux  (taken  positive  downward),  however,  varies  linearly 
with  the  SST  with  a slope  of  about  -105  cal  cm”^  day”*  °K-1 
and  accounts  for  practically  all  of  the  variation  of  the  net 
surface  heat  flux  with  respect  to  SST.  The  surface  evapora- 
tion is  found  to  increase  wit*1,  increasing  cloud  cover  at  a 
fixed  SST,  but  its  variation  with  respect  to  SST  is  inde- 
pendent of  small  changes  in  the  cloud  cover.  In  addition, 
the  cloud  cover  increases  with  SST  if  all  other  imposed 
parameters  are  the  same. 

From  the  standpoint  of  OTPP  influence  on  regional 
climate,  the  change  induced  in  the  surface  evaporation  rate 
is  the  most  significant  quantity  since  it  is  probably  well 
correlated  with  rainfall.  The  present  study  indicates  that 
the  surface  evaporation  rate  changes  substantially  with  the 
SST  and,  consequently,  must  be  given  proper  consideration 
when  planning  the  deployment  of  OTPP's.  For  example,  if 


a 
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the  maximum  tolerable  change  in  the  surface  evaporation 
rate  is  taken  to  be  -5%,  we  find  that  maximum  allowable 
SST  perturbation  is  about  -0.2°  C. 

The  present  study  can  be  extended  in  several  ways. 
First  of  all,  the  cumulus  parameterization  model  can  be 
improved  by  dropping  the  assumption  that  the  cloud  updrafts 
begin  with  the  mixed  layer  values  of  dry  static  energy  and 
water  vapor  mixing  ratio  (see  Equations  (55)  and  (56)). 
Rather,  following  Ogura  and  Cho  (1974),  we  can  recognize 
the  fact  that  updrafts  in  the  mixed  layer  which  ultimately 
produce  cumulus  clouds  tend  to  be  slightly  warmer  and  more 
moist  than  their  surroundings.  This  will  allow  the 
existence  of  cloud-induced  eddy  fluxes  of  heat  and  moisture 
out  of  the  mixed  layer.  Note  that  these  eddy  fluxes  do  not 
occur  in  the  present  model  since  the  rates  at  which  heat 
and  moisture  are  transported  out  of  the  mixed  layer  into 
rising  cumulus  clouds  are  exactly  balanced  by  the  rates  at 
which  these  quantities  are  put  back  into  the  mixed  layer  by 
compensating  subsidence  in  the  environment. 

In  addition,  the  case  of  saturated  mixed  layers 
should  be  accounted  for.  Recall  that  the  present  study 
applies  only  to  the  tupical  trade  wind  situation  in  which 
cloud  base  is  slightly  above  the  top  of  the  mixed  layer. 

This  gives  rise  to  scattered  trade  wind  cumuli  which  are 
accounted  for  by  the  Arakawa-Schubert  theory.  As  discussed 
by  Schubert  (1976),  however,  cloud  base  often  falls  below 
the  top  of  the  mixed  layer  in  subtropical  latitudes.  In 
this  case  a solid  deck  of  stratocumulus  clouds  exists  in 
the  upper  part  of  the  mixed  layer  and  must  be  accounted  for 
by  another  theory.  This  could  be  extremely  important  to  the 
OTPP  problem.  For  example,  if  OTPP  operation  causes  a 
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transition  from  the  scattered  cumuli  regime  to  the  strato- 
cumulus  regime,  the  surface  radiative  fluxes  may  be  altered 
drastically.  Unfortunately  there  is  at  present  no  adequate 
theory  for  the  transition  from  an  unsaturated  to  a saturated 
mixed  layer.  However,  a forthcoming  paper  by  Randall  and 
Arakawa  (1978)  should  be  helpful  in  this  regard. 

A more  serious  drawback  of  the  present  study  is 
that  it  treats  the  ocean  as  an  infinite  heat  reservoir  and 
ignores  the  highly  interactive  nature  of  the  upper  ocean 
and  lower  atmosphere.  A more  realistic  approach  would  be 
to  couple  the  marine  boundary  layer  model  described  here  to 
an  appropriate  upper  ocean  model.  Then,  OTPP's  could  be 
treated  as  diabatic  heat  sinks  and  the  coupled  model  could 
be  used  to  predict  the  resulting  equilibrium  SST  and  the 
surface  fluxes  for  various  degrees  of  OTPP  activity. 

A coupled  model  as  envisioned  above  would,  of 
course,  be  useful  in  air-sea  interaction  studies  other  than 
those  related  to  the  OTPP  problem.  One  possible  application 
of  such  a sophisticated  model  would  be  to  test  the  validity 
of  simpler  models  linking  the  ocean  and  atmosphere  for  the 
purpose  of  atmospheric  and/or  oceanic  forecasting. 

A coupled  air-sea  model  of  this  type  should 
eventually  be  extended  to  two  dimensions.  This  would  allow 
an  explicit  treatment  of  horizontal  advection  and  the 
characteristics  of  the  atmosphere  downstream  from  the  OTPP- 
generated  SST  anomaly  could  be  examined.  In  addition,  the 
large-scale  vertical  motion  and  one  component  of  the 
geostrophic  wind  could  be  predicted  rather  than  prescribed, 
thereby  allowing  the  SST  anomaly  to  alter  the  large-scale 
flow  pattern. 
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Appendix  A 
LIST  OF  SYMBOLS 


A(A) 


Cloud  work  function  for  subensemble  X. 


ALg(A)  Cloud  work  function  for  subensemble  X 

generated  by  large-scale  processes  alone  in 
one  time  step  of  the  model. 

a Albedo  of  air  column  from  sea  surface  to 

° model  top. 


ag  Albedo  of  sea  surface. 

B Stefan-Boltzman  constant. 


c 


Total  condensation  rate  per  unit  mass  of  air. 


Specific  heat  at  constant  pressure  for  dry 
air. 


c(  A , p ) 
D(p) 

e 

Eo 

aV^o 


f 

F(  A ) 

f£lr,  fcld 


Raindrop  generation  rate  parameter. 

Condensation  rate  for  subensemble  A at  level 
p per  unit  cloud  nass  flux. 

Total  detrainment  rate  for  the  cloud  ensemble 
at  level  p. 

Total  evaporation  rate  per  unit  mass  of  air. 

Net  heat  flux  at  the  sea  surface  (positive 
downward) . 

Change  in  the  net  surface  heat  flux  with 
respect  to  sea  surface  temperature. 

Coriolis  parameter. 

Large-scale  forcing  for  subensemble  A.  * 

Net  downward  flux  of  long-wave  (infrared) 
radiation  in  clear  region,  cloudy  region. 


F*  Net  downward  flux  of  long-wave  (infrared) 

° radiation  at  sea  surface. 


A-l 


g 


H 

h 

h* 

h0(X,p) 

Ho 

k 

K* 

K(  X , X " ) 
Ky( X , X ' ) 
Kd(X,X') 

km(X) 

Lv 

L 

i(X,p) 


Acceleration  of  gravity. 

Atmospheric  scale  height,  RT/g. 

Moist  static  energy. 

Saturated  moist  static  energy. 

Moist  static  energy  of  clouds  in  subensemble 
X at  level  p. 

Surface  sensible  heat  flux  (positive  downward) 
Entrainment  parameter 

Vertical  eddy  diffusion  coefficient  for 
momentum . 

Kernel  from  cumulus  parameterization  theory. 
Vertical  mass  flux  kernel. 

Detrainment  kernel. 

Mixed  layer  kernel. 

Latent  heat  of  evaporation  for  water. 
Monin-Obukhov  length. 

Liquid  water  mixing  ratio  of  clouds  in  sub- 
ensemble X at  level  p (mass  of  liquid  water 
per  mass  of  dry  air). 

Liquid  water  mixing  ratio  of  cloud  air 
detrained  at  level  p. 

Surface  latent  heat  flux  (positive  downward). 

Vertical  mass  flux  (positive  upward)  at  pb 
for  subensemble  X. 

Total  number  of  clouds  in  cumulus  ensemble. 

Ocean  Thermal  Power  Plant. 

Pressure,  reference  pressure. 

Pressure  at  surface,  first  gridpoint  above 
surface,  top  of  surface  layer. 


f 


( 


f 


r 


i 

‘ 

! 

€ 


« 


C 


t 


y P 


> 


1 


V V pc 


Pd(X) 

pe 

Qr 

q 

qm 


Q0(X,p) 


q* 

q(p) 


qo’  q10 


Pressure  at  top  of  mixed  layer,  top  of 
transition  layer,  cloud  base. 

Pressure  at  level  of  detrainment  (level  of 
cloud  tops)  for  clouds  in  subensemble  A. 

Pressure  at  model  top. 

Net  radiational  heating  rate. 

Water  vapor  mixing  ratio  (mass  of  water  vapor 
per  mass  of  dry  air). 

Horizontally-averaged  water  vapor  mixing  ratio 
of  mixed  layer. 

Water  vapor  mixing  ratio  of  clouds  in  sub- 
ensemble A at  level  p. 

Saturation  value  of  water  vapor  mixing  ratio. 

Water  vapor  mixing  ratio  of  cloud  air 

Water  vapor  mixing  ratio  at  sea  surface, 

10  m height . 


-Qo 

Surface  evaporation  rate. 

R 

Gas  constant  for  dry  air. 

R* 

Total  raindrop  generation 
of  air. 

rate 

per  unit  mass 

r( A ,p) 

Raindrop  generation  rate 
at  level  p per  unit  cloud 

for  subensemble  A 
mass  flux. 

s 


Dry  static  energy. 


s 


m 


Horixontally-averaged  dry  static  energy  of 
mixed  layer. 


s0( A,p) 
svo(X*p) 


Virtual  dry  static  energy. 

Dry  static  energy  of  clouds  in  subensemble  A 
at  level  p. 

Virtual  dry  static  energy  of  clouds  in  sub- 
ensemble A at  level  p. 
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C 


Dry  static  energy  of  cloud  air  detrained  at 
level  p. 


s(p) 


s(p) 


SCLR,  SCLD 


ss 


SST 

t 

T. Tj 

To 

u 

Ug 

U, 


u,  v1 


Uc 

V 

V 


g 


V 

-►og 


Incident  short-wave  (solar)  radiation  at 
model  top. 

Downward  flux  of  short-wave  radiation  subject 
to  scattering. 

Downward  flux  of  short-wave  radiation  subject 
to  absorption. 

Dry  static  energy  of  cloud  air  detrained  at 
level  p. 

Net  downward  flux  of  short-wave  (solar) 
radiation  in  clear  region,  cloudy  region. 

Net  downward  flux  of  short-wave  (solar) 
radiation  at  sea  surface. 

Sea  Surface  Temperature. 

Time. 

Temperature,  reference  temperature 

Temperature  at  sea  surface. 

Component  of  wind  in  x-direction 

Component  of  geostrophic  wind  in  x-direction. 

Surface  atmospheric  friction  velocity. 

Wind  speed,  wind  speed  at  first  gridpoint 
above  sea  surface. 

Wind  speed  at  sea  surface. 

Component  of  wind  in  y-direction. 

Component  of  geostrophic  wind  in  y-direction. 

Wind  velocity  vector  on  a surface  of  constant 
pressure . 

Geostrophic  wind  velocity  vector  at  sea 
surface . 


Fraction  of  model  domain  covered  by  i-th 
cloud. 


r 

°i 

Fraction  of  model  domain  covered  by  i-th 
cloud. 

K 

von  Karmen's  constant. 

X 

Entrainment  rate  for  clouds  in  subensemble 

r 

xd(p) 

Entrainment  rate  for  clouds  that  have  level 
of  detrainment  (level  of  cloud  tops)  at 
level  p. 

r 

n( X ,p) 

Ratio  of  vertical  mass  flux  due  to  subensemble 
X at  level  p to  vertical  mass  flux  due  to 
subensemble  X at  pb- 

to 

Vertical  p-velocity,  dp/dt. 

c 

fa> 

a 

Average  vertical  p-velocity  in  cloud  ensemble 
times  fraction  of  model  domain  covered  by  the 
ensemble . 

el 

Q 

>* 

tJ 

Average  vertical  p-velocity  in  type  X clouds 
at  level  p times  fraction  of  model  domain 
covered  by  subensemble 

It 

u> 

0 b 

Average  vertical  p-velocity  in  cloud  ensemble 
at  pb  times  fraction  of  model  domain  covered 
by  the  ensemble. 

“b 

Vertical  p-velocity  at  pb  averaged  across 
entire  model  domain. 

< 

*b 

Vertical  p-velocity  at  pb  averaged  across  part 
of  model  domain  not  occupied  by  clouds. 

To 

Surface  wind  stress  vector. 

If 

T*  Ty 
o o 

Component  of  surface  wind  stress  vector  in  x- 
direction,  y-direction. 

At 

Time  step  in  numerical  integration. 

L 

As 

Jump  of  dry  static  energy  across  transition 
layer . 

Aq 

Jump  of  water  vapor  mixing  ratio  across 
transition  layer. 

k 

A8v 

Jump  of  virtual  dry  static  energy  across 
transition  layer. 

if 


A-6 


Unit  vector  in  x-direction,  y-direction. 

Gradient  operator  at  constant  pressure. 

Horizontal  average  on  constant  pressure 
surface  taken  across  entire  model  domain. 

Horizontal  average  on  constant  pressure 
surface  taken  across  part  of  model  domain 
not  occupied  by  clouds. 

Deviation  from  horizontal  average  taken  across 
entire  model  domain. 

Average  taken  over  depth  of  mixed  layer. 

Vertical  eddy  flux  of  dry  static  energy  at 
surface,  pb . 


Vertical  eddy  flux  water  vapor  mixing  ratio 
at  surface,  pfa . 


Vertical  eddy  flux  cf  virtual  dry  static 
energy  at  surface,  p^. 


Total  vertical  cloud-induced  flux  of  dry 
static  energy,  water  vapor  mixing  ratio, 
liquid  water  mixing  ratio. 
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Appendix  B 


DETERMINATION  OF  SATURATION  MIXING 
RATIO  AND  GAMMA 


By  definition,  the  saturation  mixing  ratio,  q*, 
is  given  by 


0.622  e (T) 

«*(T-P)  ' P - ».*T) 


(B.l) 


s 


where  e (T)  is  the  saturation  vapor  pressure  over  a plane 
s 


surface  of  water.  Following  Washington  and  Kasahara 
( 1970) , we  take 


e (T)  = 6.11  exp 

S 


17. 269(T  - 273.16) 


T - 35.86 


(B.  2) 


where  e (T)  is  in  mb  and  T is  in  degrees  Kelvin, 
s' 


The  function  y is  defined  by 


1 


Appendix  C 

SOLUTION  OF  THE 
FREDHOLM  INTEGRAL  EQUATION 

The  discrete  representation  of  (65)  is 
m 

G‘V  - I 

0*1 

for  i * 1 ,2 , . . . , n , (C.l) 


Ki,J  QJ  Ak 


where  depends  on  the  quadrature  formula  used. 

AX  = X.+-  - X.,  and  n is  the  number  of  cloud  types  repre- 
0 * 0 

sented  in  the  model.  An  iterative  scheme  for  solving 
(C.l)  is 


- °(h) 


where  A is  an  acceleration  factor.  Equation  (C.2)  can  be 
written  in  fully  implicit  form  as  a matrix  equation 


Equation  (C.3)  can  be  rearranged  to  give 


m 


(k+1) 


ALS  D (k)  m(k) 
At"  + 2 m + ~ 


where 


and  I is  the  identity  matrix.  Equation  (C.4)  is  solved 
iteratively  for  the  elements  of  m which  yield  the  best 
solution  of  (C.l)  subject  to  the  constraints  of  (65). 
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